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Abstract 

We show that the previously introduced CKMT model, based on conventional Regge 
theory, gives a good description of the HERA data on the structure function for large 
rapidity gap (diffractive) events. These data allow, not only to determine the valence and 
sea quark content of the Pomeron, but also, through their Q 2 dependence, give information 
on its gluonic content. Using DGLAP evolution, we find that the gluon distribution in 
the Pomeron is very hard and the gluons carry more momentum than the quarks. This 
indicates that the Pomeron, unlike ordinary hadrons, is a mostly gluonic object. With our 
definition of the Pomeron flux factor the total momentum carried by quarks and gluons 
turns out to be 0.3 0.4 - strongly violating the momentum sum rule. 
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1. Introduction 

Diffraction dissociation of virtual photons, observed by the HI and ZEUS collabora- 
tions at HERA [1-4] provides information on the nature of the Pomeron and on its partonic 
structure. In a previous paper [5] a model (CKMT model) based on Regge theory has been 
proposed for the description of diffractive dissociation of both real and virtual photons. 
In this model the structure function of the Pomeron Fp((3,Q 2 ) is related to that of the 
deuteron [6] via Regge factorization. Recently, experimental information on the Pomeron 
structure function at different values of Q 2 has been obtained [2, 4]. These data are in 
reasonable agreement with the predictions of the model. 

In this paper, we extend our analysis of the Pomeron structure function and determine 
the gluonic content of the Pomeron. We use our model as initial condition in DGLAP 
evolution. We show that the Q 2 -dependence of the diffractive structure function, , 
observed experimentally, can only be understood if the gluon distribution in the Pomeron 
is very hard. It is also shown that the gluons carry more momentum than quarks. This 
observation confirms theoretical expectations that the Pomeron is mainly connected with 
gluonic degrees of freedom. 

The plan of the paper is as follows. In Section 2 we recall the ingredients of the 
CKMT model. In Section 3 we perform the QCD evolution and in Section 4 we present 
our numerical results. In Section 5 we discuss other recent QCD analysis and give our 
conclusions. 

2. - The CKMT model 

The contribution of diffractive processes to the deep inelastic scattering (DIS) cross- 
section, corresponding to the Pomeron exchange diagram of Figs. 1 and 2, can be written 
in the form [4] 




where x and y are standard DIS variables, t is the invariant momentum transfer t = 



— (p — p 1 ) 2 and xp is the fraction of the proton momentum carried by the Pomeron : 

= VP - Vl M 2 + Q 2 = x_ 
XP Pp-Pj ~ W 2 + Q 2 ~ (3 ' U 

We have argued in Ref. [5] that the Pomeron in the diagrams of Figs. 1 and 2 can be con- 
sidered as a Regge pole with a trajectory a pit) = ap(0) -\-a't determined from an analysis 
of soft processes, in which absorptive correction (Regge cuts) are taken into account. We 
have [5], [7] 

a P (0) = 1.13 , a'p(O) = 0.25 GeV" 2 . (3) 
In this case the diffractive contribution to DIS can be written in a factorized form 

, 2 
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where g pp it) is the Pomeron-proton coupling g pp {t) = g p (0)exp(Ct) with {g pp {0)) 2 = 23 
mb and C = 2.2 GeV -2 . Eq. (4) is the definition* of the Pomeron structure function 
Fp(f3,Q 2 ,t) with the variable (3 = x/xp playing the same role as the Bjorken variable x 
in F% . As emphasized in Ref. [5], the definition of Fp depends on the particular choice 
of the Pomeron flux factor. For instance in Ref. [8] the flux factor differs from ours by a 
factor 2/tt. (See also Ref. [9]). 
At large Q 2 we have 

F P (p, Q 2 ,t) = J2 e 2 iP [qf{P, Q 2 ,*) + qf{P, Q 2 ,*)] ■ (5) 

i 

Due to the arbitrariness in the normalization of Fp discussed above, the partonic distri- 
butions qf and qf do not satisfy, in general, the energy-momentum sum rule. 

Many models of the Pomeron structure function are available in the literature [8,10- 
14]. In the CKMT model [5] Fp is determined using Regge factorization together with the 



* Note that the definition of Fp in Eq. (4) differs from the one in Ref. [5] by a factor 
(1 — (3). The definition in Eq. (4), which coincides with the one used in experimental 
papers [2, 4], is more appropriate for a partonic interpretation. We thank X. Artru for an 
enlightening correspondence on this point. 
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values of the triple Regge couplings determined from soft diffraction data. More precisely, 
the proton and Pomeron structure functions are given in terms of the diagrams of Figs. 3 
and 4, respectively. (Note that Fig. 4 is the upper part of the diagram obtained by squaring, 
in the sense of unitarity, the 7P amplitude in Fig. 1). We see that, due to factorization, 
F2 can be obtained from (or more precisely from the combination F$ = \{F% + F£)), 
by replacing the Reggeon-proton couplings and gf appearing in the lower part of 
Fig. 3, by the corresponding triple Reggeon couplings rppp and rppf which appear in the 
lower part of Fig. 4. The latter can be determined from soft diffraction data (see below). 
Another difference between the two structure functions is the (3 — > 1 behaviour. Indeed, 
as far as the dimensional counting rules are concerned, the Pomeron is like a pion and, 
therefore, there is one spectator less than in the proton case. Thus the power of 1 — (3 is 
smaller by two units in the Pomeron case. (The same result is obtained at Q 2 = in terms 
of Regge intercepts). The CKMT model can be summarized with the following formulae 
[5, 6], valid in the region 1 < Q 2 < 5 GeV 2 : 

Fp((3,Q 2 ) = F* ((3,Q 2 ;A^ eA,B fB,n n - 2) , (6) 

F 2 d (x, Q 2 ) = A(Q 2 ) x~ A( ^ (1 - x) n( ^+ 4 + B x 1 -^ (1 - x) n( ^ (7) 

where e = fpp(0)/gp p (0) and / = rpp(0)/gf (0). The values of all other parameters 
are given in Refs. [5] [6]. A main feature of the CKMT model is the Q 2 dependence of 
the effective Pomeron intercept, A = a e J^ — 1, which appears in the upper part of the 
diagrams in Figs. 3 and 4. It was argued in Ref. [6] that this is due to the fact that the 
size of the absorptive corrections decreases when Q 2 increases. As a consequence, when 
using a parametrization with a simple power x _A , A must depend on Q 2 . 

The parameters e and / can be determined from soft diffraction data. Here also 
absorptive corrections, which are very important in diffractive processes [15], have to be 
taken into account when extracting the values of e and / from the data. In an analysis 
of soft diffraction without absorption [16] one obtains e = f = 0.025. In the first paper of 
Ref. [15] it was shown that these values have to be multiplied by a factor of about three in 
order to take into account absorptive effects. The values of e and f have, of course, some 







uncertainties-the largest being in the value of / [5]. Experimentally, the t-dependence of 
the triple Reggeon couplings is very small. We have taken it to be the same for r-p P and 
rppf and included it in the function C. In this way the t-dependence of Fp, which is 
expected to be very small, is factored out. 

The comparison of the prediction [5] of the CKMT model with recent data from the HI 
and ZEUS collaborations is quite satisfactory. The data confirm the factorized dependence 
of the hard diffractive cross-section on xp (Eq. (1)) with the parameter of the Pomeron 
trajectory in Eq. (2). Agreement in absolute values is also good for most values of (3 and 
Q 2 , favoring our hypothesis that soft and hard diffraction are governed by the same triple 
Regge couplings. 

3. QCD evolution 

A closer look at the model predictions shows however a systematic decrease with Q 2 
at (3 > 0.2. Such a trend is not present in the data. The reason for such a decrease is the 
following. As explained in Ref. [5] the model predictions at large Q 2 have been obtained 
from DGLAP evolution [17], using Eqs. (6) and (7) at Qq as initial condition. The result 
depends, of course, on the input gluon distribution function. In Ref. [5] it was assumed 
that the relation between gluon and sea quark distributions in the Pomeron was the same 
as in the proton. This corresponds to a rather soft gluon distribution in (1 — (3) 3 at Q 2 
5 GeV 2 . It was already pointed out in Ref. [5] that this assumption was not justified and 
should be changed. Indeed, such a behaviour is valid for ordinary mesons with a large 
valence quark content. However, the Pomeron can have a much harder gluon distribution 
("valence gluons"). This is also in agreement with UA8 data [18]. 

In the following we will perform DGLAP Q 2 -evolution using the valence and sea 
quark initial distribution given by Eqs. (6) and (7) at Qq = 5 GeV 2 , with a harder gluon 
distribution. More precisely, we use the same (3 — > behaviour as for sea quark [5], but 
we leave the power n g of 1 — (3 as a free parameter, i.e. we put 



where A g = 1.71 is the gluon normalization constant in a nucleon [6] and e and A are 
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defined in Eq. (7). Note that the gluon normalization is obtained from Regge factorization 
and the energy-momentum sum rule is not used. For n g small enough, the Q 2 behaviour 
of Fp will be very different from the one of F^- The latter exhibits a decrease with Q 2 for 
x > 0.2 due to a softening of the quark distribution resulting from gluon emission. Even 
though the same effect is present for the Pomeron, a harder distribution of gluons will lead 
in this case, through gluon decay, to an increase with Q 2 in the number of quarks up to 
comparatively large values of (3. 

Numerical calculations have been performed using the QCD evolution program of 
Refs. [19, 20, 14]. We present the results in a one loop approximaton but we have checked 
that practically the same results are obtained in two loops. The input and evolved partonic 
distributions are given in Figs. 5-8. 

4. Numerical results 

We have found [21] that a good description of the Q 2 dependence of the HERA data 
is achieved with either n g = or — 0.5 at Qq = 5 GeV 2 . We have also allowed the 
parameters e and / in Eqs. (6) and (8) to vary within limits consistent with the results 
from the analysis of soft diffraction. A good description of the data is obtained with 
e = f = 0.07. All other parameters in Eq. (7) have the values given in Ref. [6]. 

The 3-dimensional structure function F^(x, Q 2 , xp), integrated over t, is plotted ver- 
sus xp in Fig. 9 and compared with HI data. The results for the Q 2 and (3 dependence of 
F® (integrated over t and xp) are given in Fig. 10 and compared with HI data. In both 
cases, unpublished data from the ZEUS collaboration are also available. These data are in 
agreement, within errors, with the ones of HI. We see that the Q 2 -dependence at large (3 
is flat, in agreement with experiment - whereas in Ref. [5] there was a substantial decrease 
with increasing Q 2 . 

In Fig. 11 we give the Pomeron structure function defined in Eq. 4 at different values 
of Q 2 . We see that the cross-over, which is at (3 ~ 0.2 when the gluon distribution in the 
Pomeron is taken to be proportional to that of the proton [5], is shifted to larger values of 
(3. Comparison of Figs. 11a (n g = 0) and lib (n g = —0.5) shows that this shift is larger 



when the gluon distribution is harder. 

In all figures, except Fig. 11, we have presented the results only for n g = —0.5. 
However, for n g = there is very little change in our results - at least in the range of 
variables where data are available. 

Finally, we give in Fig. 12 the fractions of the Pomeron momentum carried by quarks 
and gluons. With n g = — 0.5 (n g = 0). For Qq = 5 GeV 2 , we have 



We see that for both values of n g the momentum fraction carried by gluons is larger than 
the one carried by quarks, indicating that the Pomeron has a substantial gluonic component 
with a very hard gluon distribution. The energy-momentum sum rule is strongly violated. 

5. Comparison with other QCD analysis and conclusions 

Very recently, QCD-based analysis of the Pomeron structure function and its Q 2 
evolution have been carried out [13, 22, 23, 24]. In these papers the initial condition 
in DGLAP evolution equation and, in particular, the gluon momentum distribution is 
different from ours. However, in all cases conclusions similar to ours are reached, namely, 
the Q 2 dependence of the data can be reproduced only if the gluon momentum distribution 
in the Pomeron is hard - the gluons carrying more momentum than the quarks. Note that 
in Refs. [22, 23] the gluon normalization relative to that of the quarks is determined from 
the momentum sum rule. This leads to a large normalization constant for gluons. It is 
known from UA8 results [18] that the normalization of partons satisfying the momentum 
sum rule leads to absolute predictions which exceed experiment by a factor close to five 
(see also Ref. [9]). 

In conclusion, recent HERA data confirm the predictions of the CKMT model of hard 
diffraction based on Regge theory, factorization, and the assumption that not only the 
Pomeron intercept but also the triple Reggeon couplings are the same in soft and hard 




(0.15) 
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diffraction. Although a QCD analysis of the Q 2 -dependence of the present data does not 
allow to extract the shape of the gluon distribution in the Pomeron, it requires a hard 
gluon distribution - the gluons carrying substantially more momentum than the quarks. 
This conclusion can be reached without reference to the momentum sum rule, using a 
normalization constant of the gluon distribution determined from Regge factorization. 
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Figure Captions 

Fig. 1 Pomeron exchange diagram for single diffraction dissociation of a virtual photon. 

Fig. 2 Pomeron exchange diagram for a double diffraction dissociation in which both the 
hard current and the target proton are excited. 

Fig. 3 Regge exchange diagram for the proton structure function. 
Fig. 4 Regge exchange diagram for the Pomeron structure function. 

Fig. 5 Initial gluon distribution (Q 2 = 5 GeV 2 ) and its perturbative QCD evolution for 
n g = —0.5. 

Fig. 6 Same as Fig. 5 for a u quark. 
Fig. 7 Same as Fig. 5 for an s-quark. 

Fig. 8 Same as Fig. 5 for a c-quark. The charm distribution is taken to be zero at Qq and 
is generated, at larger Q 2 , via DGLAP evolution. 

Fig. 9 The 3-dimensional proton structure function for diffractive events F^ '(x,xp ,Q 2 ) is 
plotted versus xp at various values of Q 2 and (3 = x/xp and compared with HI data 
[4]. The model results are for n g = —0.5. Very similar results are obtained with 
n g = 0. 

Fig. 10 The 2-dimensional proton structure function for diffractive events F^(J3, Q 2 ) is plot- 
ted versus Q 2 at fixed (3 and versus (3 at fixed Q 2 and compared with HI data [4]. The 
model results are for n g = —0.5. Very similar results are obtained with n g = 0. In or- 
der to take into account the presence of double diffraction in the data, the theoretical 
result has been multiplied by 1.3. 

Fig. 11 The Pomeron structure function F^"(J3, Q 2 ), defined in Eq. (4), is plotted versus (3 at 
several values of Q 2 for n g = —0.5 (Fig. 11a) and n g = (Fig. lib). 

Fig. 12 The functions of the Pomeron momentum carried by quarks (N q ) and gluons (N g ) 
defined by eqs. (9) and (10) are plotted versus Q 2 for two different values of n g . 
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